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Abstract. M0S2 and WS2 layered transition-metal dichalcogenides are indirect band gap semiconductors 
in their bulk forms. Thinned to a monolayer, they undergo a transition and become direct band gap 
materials. Layered structures of that kind can be folded to form nanotubes. We present here the electronic 
structure comparison between bulk, monolayered and tubular forms of transition metal disulfides using 
first-principle calculations. Our results show that armchair nanotubes remain indirect gap semiconductors, 
similar to the bulk system, while the zigzag nanotubes, like a monolayer, are direct gap materials, what 
suggests interesting potential applications in optoelectronics. 

PACS. PACS-key discribing text of that key - PACS-key discribing text of that key 



^ 1 Introduction 



The discovery of carbon nanotubes^ immediately stim- 
ulated the research to explore other inorganic materials, 
which may form tubular structures. Nanoparticles of var- 
ious inorganic layered materials, such as WS2 or M0S2, in 
analogy to carbon, can form nanostructures of fullerene- 
like and nanotubular shapes. In 1992 and 1993, Tenne 
and co-workers have shown that layered transition-metal 
dichalcogenides (LTMDCs), namely WS2 and M0S2, form 
the so-called inorganic nanotubes and fullerene-like nano- 
particles. Nowadays, it is well-known that several lay- 
ered inorganic compounds posses structures comparable 
to graphite (honey-comb-like) and as such can form tubes. 
Among them are transition metal dichalcogenides, halides, 
and oxides. 

LTMDCs of TMX2 type (TM = Mo, W, Nb, Re, Ti, 
etc., X = S, Se, Te) have been extensively studied ex- 
perimentally and theoretically for the last 40 years.—"— 
Molybdenum disulfide is a prototypical LTMDC, which 
is composed of two-dimensional trilayered S-Mo-S sheets 
stacked on top of one another and held together by weak 
bonds (see Fig. [T] left). Each sheet is trilayered with a 
Mo atom in the middle that is covalently bonded to six 
S atoms located in the top and bottom of the sheet (see 
Fig.HD. 



LTMDCs are of special interest because of their ex- 
traordinary properties and diverse applications. We have 
recently shown that TMS2 (TM = Mo, W) are indirect 
band gap materials in their bulk form and they become 
direct band gap semiconductors when thinned to mono- 
layers.— These results support the recent experimental 
findings for M0S2— and show that quantum confinement 
is the key factor to tune electronic structures of LTMDCs. 
Therefore, electronic and optical properties of LTMDCs 
are very interesting for fabrication of field-effect transis- 
tors, as it was recently realized by Kis and co-workers 
for a M0S2 monolayer, -'^-d''' In this transistor, Hf02 was 
used as a gate insulator.— It exhibits a room-temperature 
current on/off ratio exceeding 1x10^ and mobility com- 
parable to the mobility of thin silicon films or graphene 
nanoribbons. Other applications, such as catalysis, opto- 
electronics and photovoltaics, have been proposed and in- 
vestigated as welL^iiiS^i^a 



Bulk LTMDCs can be synthesized using a variety of 
methods. Generally, dichalcogenides are prepared by heat- 
ing under hydrogen sulfide flow. Monolayers of such lay- 
ered materials can be produced by e.g. liquid exfoliation, 
as successfully performed by Coleman et al.— Moreover, 
it was shown that it is possible to distinguish the number 
of layers of dichalcogenide crystals, using a simple optical 
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Fig. 1. (Color online) The atomic structures of layered (left) and tubular (middle - zigzag and right - armchair) transition-metal 
disulfides of TMS2 type (TM = Mo, W). (top) the cross sections and (bottom) the hexagonal sheets are shown for both types 
of the systems. 



model, by calculating the contrast of nanolayers deposited 
on Si02 wafer with varying thickness. 

Transition metal dichalcogenide nanotubes (TMDC- 
NTs) can be produced using e.g. chemical vapor transport 
technique^ or by high-temperature annealing of the re- 
spective metal trisulfides.^^ For more details, refer to the 
comprehensive review on inorganic nanotubes by Nath et 

TMDC-NTs have been a subject of various investiga- 
tions^i2Sr— in order to understand physical and chemical 
properties of these materials and to propose important ap- 
plications. For example, it has been shown that TMDCs- 
NTs behave as exceptional lubricants. ^^Sj^ The mechani- 
cal properties of WS2 nanotubes under axial tension and 
compression,'^^ and M0S2 nanotubes under squeezing,— 
have been investigated in terms of the lubrication pro- 
cess at high loads. Experimental results show that at high 
loads the lubricant is expelled from the location of close 
interface distance due to mechanical pressure, and it hence 
looses its tribological properties, leading to an increased 
friction coefficient. If the M0S2 nanotubes or -onions are 
added to base grease, the friction coefficient remains low, 
even at very high loads. Density functional-based calcu- 
lations show that under the squeezing of nanotubes the 



M0S2 platelets are formed, partially attached to the grips, 
which provide good lubrication at the position of closest 
contact of the materials sliding against each other. This 
excellent lubrication of nanostructures is interpreted as 
'nano-coating'. The mechanical behaviour of WS2 nan- 
otubes under axial tension and compression— shows that 
they are ultra-strong and elastic, what distinguishes them 
from other known materials. The failure of the nanotubes 
is abrupt starting at a single atomic defect and propagat- 
ing very quickly across its entire circumference. Moreover, 
M0S2-NTS have been recently used for catalytic conver- 
sion of carbon oxide and hydrogen into methane and wa- 
ter. These foundings are quite unexpected, as the fully 
bonded sulfur atoms in the LTMDC surfaces are not ex- 
pected to be chemically active. 

Doping these semiconducting nanotubes may lead to 
new optoelectronic nanomaterials. Ivanovskaya et al.— 
have investigated the effect of Mo to Nb substitution on 
the electronic structure of M0S2 nanotubes using density 
functional based tight-binding (DFTB) method.^i^ It 
has been found that composite Moi_a;Nba;S2 nanotubes 
(with Nb contents of 5, 10 and 25 at%) are more stable 
than the corresponding pure tubes. This effect was even 
stronger for larger tube diameters. At room temperature 



Nourdine Zibouche et al.: From Layers to Nanotubes: Transition Metal Disulfides TMS2 



3 



the Nb dopant prefers to arrange in ordered manner along 
the tube axis, whereas the site entropy contribution favors 
a random distribution at high temperatures. Ivanovskaya 
et al.— reportred that all the doped nanotubes studied in 
this work have metallic properties, independent of their 
chirality, diameters or the substitutional patterns. The 
density of states (DOS) close to the Fermi level of Nb- 
substituted M0S2 nanotubes can be tuned in a wide range 
by the degree of doping. On the other hand, Ivanovskaya 
et al.— have also studied the Nb substitution in M0S2 
layered materials and found out that the Nb atoms prefer 
to distribute homogeneously. 

The properties of inorganic nanotubes can be further 
widely explored, as they are much less investigated as 
their carbon counterparts. Especially, understanding of 
the change in electronic properties, when going from lay- 
ers to nanotubes or after substitution, is of high interest 
due to potential applications in nano- and optoelectronics. 

In this paper, we would like to present the compari- 
son of the electronic structure between bulk, monolayered 
and tubular forms of TMS2, where TM = Mo and W (see 
Fig. [T] middle and right). The first principle calculations 
were performed using localized Gaussian basis functions 
and compared to the available experimental and theo- 
retical data. Our results show that TMDC-NTs of arm- 
chair type are indirect band gap semiconductors (like a 
bulk structure), while their zigzag counterparts have di- 
rect band gaps and resemble the electronic structure of a 
TMDC monolayer. 

2 Methods 

In this work, we have studied TMDC-NTs of TMS2 
type (where TM = Mo and W) with different tube sizes for 
both armchair and zigzag configurations. All layered struc- 
tures, initial systems to built up nanotubes, have hexag- 
onal symmetry and belong to the P63 / mmc space group. 
The monolayers were cut out from the fully optimized bulk 
structures as (0 1) surfaces. 

First-principle calculations were performed on the ba- 
sis of density functional theory (DFT) as implemented 
in the CRYSTAL09 code.— The exchange and correla- 
tion terms were described using general gradient approx- 
imation (GGA) in the scheme of PBE (Perdew-Burke- 
Ernzerhof). ■^^ The following basis sets were used: Mo_SC_- 
HAYWSC-311(d31)G.cora_1997 (for Mo atoms),^ W_co- 
ra.l996 (for W atoms), and S_86-311G*Jichanot.l993 
(for S atoms).**- We have already shown, that the com- 
bination of these bases and PBE functional gives very 



Table 1. Calculated and experimental lattice parameters [A] 
of 2D hexagonal transition-metal dichalcogenides in the form 
of TMS2 (TM = Mo, W). Results obtained at the DFT/PBE 
level. 



Structure 


Theory 


Exp. 


4,5,43 


a 


c 


a 


c 


M0S2 


3.173 


12.696 


3.160 


12.295 


WS2 


3.164 


12.473 


3.154 


12.362 



good results for electronic structure calculations of lay- 
ered TMDCs.i5 

The shrinking factor for bulk and layered structures 
was set to 8, what results in the corresponding number 
of 50 and 30 /c-points in the irreducible Brillouin zone, 
respectively. The tubular structures were treated with the 
shrinking factor 8 and 4, that correspond to 5 and 3 k- 
points along the tube axis, respectively. The mesh of k- 
points was obtained according to the scheme proposed by 
Monkhorst and Pack.— Band structures were calculated 
along the high symmetry points using the F — M — K — F 
and F — K — F paths for bulk/monolayer and nanotubes, 
respectively. 

The bulk, layers and nanotube calculations have been 
carried out employing 3D, 2D and ID periodic boundary 
conditions, respectively. 

Optimization of initial experimental structures was per- 
formed using analytical energy gradients with respect to 
atomic coordinates and unit cell parameters within a quasi- 
Newton scheme combined with the BFGS (Broyden-Flet- 
cher-Goldfarb-Shanno) scheme for Hessian updating. The 
optimized lattice parameters for all the studied materials 
are given in Table [TJ 

3 Results and Discussion 
3.1 Structural Properties 

Our calculations were performed for prototypical (n,0) 
zigzag and (n,n) armchair TMS2 nanotubes (TM = Mo 
and W) as function of n. The tube diameters {d) range 
between 12.4 A and 24.9 A for zigzag NTs, and 19.2 A 
and 42.3 A for armchair ones. These correspond to the 
index n = 11-24. 

We have compared the structural properties of dif- 
ferent NTs with those of bulk/mo nolayered systems (see 
Fig. [5]) . The lattice parameter a changes for zigzag forms 
and we found an increase by around 0.14 A when increas- 
ing d. Generally, we find the Mo-S bonds longer that their 
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W-S analogues (up to 0.1 A depending on the bond type). 
The zigzag topologies provide longer bonds (up to 0.3 A 
depending on the tube diameter) than the armchair ones. 
Layered structures have one type of TM-S bond lengths, 
while in the tubular form we can distinguish the inner 
and outer TM-S bonds (hereafter, we will refer to them as 
TM-Si and TM-So, respectively, unless otherwise stated). 



MoSg {n,n) 



MoSg {n,0) 
WSg {n,0) A-^ WSg (n,n 

3.6 



MOg layers 
WSg layers 
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Fig. 2. (Color online) Selected bond lengths of zigzag and 
armchair TMS2-NTS (TM = Mo, W) as function of the chiral 
index n. 



In detail: The optimized bond lengths of tubes are 
slightly different than those of the planar sheets. Rolling 
up a monolayer into a nanotube causes elongation of the 
TM-So bond lengths by around 0.05-0.14 A for zigzag 
and 0.01-0.02 A for armchair types (see Tab. H]). On the 
other hand, the TM-S^ bonds shrink by 0.01 A for both 
zigzag and armchair forms. TM-So bonds are longer than 



Table 2. Calculated TM-S bond lengths [A] of nanotubular 
and 2D hexagonal TMS2 (TM = Mo, W). Only the largest and 
the smallest diameter tubes are considered. For more details 
see text and Fig. [J] 



Bond 


(11,0) 


(24,0) 


(11,11) 


(24,24) 


2D 


Mo-So 


2.556 


2.466 


2.439 


2.429 


2.417 


Mo-S, 


2.402 


2.395 


2.388 


2.404 


2.417 


W-So 


2.539 


2.452 


2.427 


2.412 


2.400 


W-S, 


2.388 


2.378 


2.360 


2.381 


2.400 



TM-S, by 0.02-0.07 A for armchair NTs and by 0.07- 
0.15 A for zigzag tubes for both disulfides (depending on 
the tube diameter). Comparing various tube sizes, these 
bond lengths behave differently. Increasing the tube di- 
ameter, the TM-Si bond lengths decrease or increase for 
armchair and zigzag forms, respectively. The TM-So bond 
lengths in zigzag TMS2-NTS change for different chirali- 
ties, namely they decrease by around 0.09 A with increas- 
ing the tube diameter from (11,0) to (24,0) NT, while they 
stay almost unchanged in armchair structures. 

The TM-TM distances (within the tube circumfer- 
ence) are generally longer for tubular forms and decrease 
with increasing the diameter of the NTs. Going from (11,0) 
to (24,0) NT, these values decrease by around 0.06 A in 
the armchair configuration and by around 0.30 A for the 
zigzag. Interestingly, the TM-TM distances (in armchair 
NTs) along the axial direction do not change with the 
tube diameter and their values are similar to those in the 
layered systems. For small tube diameters, the S-S bonds 
are shorter than the corresponding bonds in the layered 
structures, but they increase with increasing d, almost ap- 
proaching the values for a bulk and a monolayer. The same 
results are found for both disulfides. 



3.2 Energetic Properties 

The calculated strain energies, i.e. the differences be- 
tween the total energies (per atom) of the tubes and the 
monolayers, scale as where d is the tube diameter 

(see Fig. [3] and Tab. [31). The correlation coefficients of the 
y = C/d'^ curves, where C is a constant, are larger than 
0.998 for all the studied tubes. The strain energy scaling is 
similar to carbon nanotubes (CNTs),— however, the en- 
ergy values are around 1 order of magnitude larger than 
for CNTs with similar diameters. We can understand this 
result easily, as it is much easier to fold a monoatomic 
layer (CNTs) than a triatomic (TMDCs). 
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Fig. 3. (Color online) Calculated strain energies of zigzag and 
armchair TMS2-NTS (TM = Mo, W) as function of the tube 
diameter (d). The energies are given per atom and related to 
the infinite monolayers of the corresponding disulfides. 
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Table 3. Correlation coefficients (r^) of the strain energy fit- 
ting as function oi C/(P for all studied tubular systems. Con- 
stant C is given in [eV A'^] and d denotes tube diameter [A]. 



NT 


C 


r^ 


Zigzag M0S2 


57.504 


0.9993 


Zigzag WS2 


58.142 


0.9990 


Armchair M0S2 


50.897 


0.9997 


Armchair WS2 


59.681 


0.9985 



The strain energies of armchair nanotubes are just 
slightly more favorable than the zigzag for a given diam- 
eter for both disulfides. This is in a good agreement with 
the work of Seifert et al.,— where the authors have used 
the DFTB method. The strain energy factor C (in eV A^) 
for armchair and zigzag M0S2 NTs on the DFTB level 
was 57.613 and 61.5206, respectively (cf. Tab. Also, 
the molybdenum forms are more stable than the tungsten 
for diameters larger than 15 A and this difference grows 
with the size of the tube. 

3.3 Electronic Properties 

We have studied the Mulliken charges [q) and the elec- 
tronic band structures of TMS2 (TM = Mo and W) in 
both layered and tubular forms. The charge transfer from 
Mo atoms to S atoms within a monolayer {qmo = 0.95, 
qs = —0.48) is similar to that in the bulk form of M0S2 
{qMo = 0.99, qs = -0.50; see Fig. ^. For M0S2, Seifert 
et al.— have found qMo = 0.90 and qs = -0.44. Slight 



-0.6 
-0.8 



' 1 ' 1 ' 1 

^.^^ 




TM-S„ _ 






1,1,1 



10 20 30 



d [A] 

Fig. 4. (Color online) Calculated Mulliken atomic charges (g) 
of tubular and layered forms of TMS2 (TM = Mo and W) as 
a function of tube diameter (d). 



differences are found for the WS2 layers, where for the 
TM atom we have the charge of 1.81 and 1.69 for the bulk 
and monolayer, respectively. The charge of sulfur atoms, 
—0.90 in bulk and —0.85 in the monolayer, is larger than 
in M0S2. 

Folding a monolayer of TMS2 into a nanotube affects 
the atomic charges mostly for small diameter sizes. Gen- 
erally, the zigzag and armchair forms have the same val- 
ues of q for a given tube diameter. The atomic charges 
approach a constant value (independent of d) for diame- 
ters larger than 25 A. Moreover, the tungsten nanotubes 
deviate from the layered structures much more than the 
molybdenum ones. 

In detail, in TMS2-NTS we should distinguish two dif- 
ferent types of sulfur atoms, S; and So, that are forming 
the inner and outer walls of the nanotubes. In case of 
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M0S2-NTS, the So atoms are slightly more negative than 
the Si atoms for both zigzag and armchair forms. Also, 
So atoms approach the atomic charges of layered struc- 
tures for larger diameters, while in the case of Si atoms, 
the q values are close to the layered structures for small- 
est tube diameters. The situation changes for WS2-NTS, 
where the So atoms are less negative than the inner ones 
and the sulfur atoms of the layered forms. The inner sul- 
furs, for WS2-NTS d smaller than 15 A, have charge close 
to the one of a bulk structure, while for d larger than 25 A 
this values are close to the charges of a monolayer. 

The charge of TM atoms becomes less positive when 
the tube diameter increases. For the smallest tube sizes, 
these values are closer to the bulk form and for d larger 
than 15 A, they approach the charges of the monolayers. 

Recently, we have studied the electronic structure of 
layered TMS2 (with TM = Mo, W, Nb, and Re).i^ Fig. [S] 
shows the band structures of M0S2 and WS2 bulk and 
monolayers. When a bulk TMS2 is thinned to a monolayer, 
there is a transition from an indirect (Z\ occurs between 
'k=r and k=l/2(/^-K)) to a direct gap semiconductor [A 
occurs at k=K). This suggests that the quantum confine- 
ment plays the key role in tailoring the electronic structure 
of such systems. 



bulk monolayer 




r MK rr mk r 



Fig. 5. (Color online) Band structures of bulk and monolay- 
ered TMS2 (TM = Mo and W). The arrows indicate the fun- 
damental band gap (direct or indirect; A) for a given system. 
The top of valence band (blue) and bottom of conduction band 
(green) are highlighted (Online color). 

Fig. [S] shows exemplary band structures around the 
Fermi level of TMS2-NTS in zigzag and armchair configu- 
rations and chiral index n = 11, 15, and 24. Band struc- 
tures of the zigzag TMS2-NTS resemble that of the cor- 



responding single- walled carbon nanotubes (SWCNTs). 
However, in contrast to SWCNTs, all the TMS2-NTS are 
semiconducting. While the band structures of (n,0) TMS2- 
NTs can be related to that of a monolayer, with direct Z\, 
the band structure of (n,n) TMS2-NTS are more related 
to the bulk, especially for larger tube diameters. 

For a given tube diameter, the band gap values of 
zigzag NTs are larger than those of armchair NTs (see 
Fig. [7]). These values, in contrast to SWCNTs, increase 
with increasing tube diameter, going from values close to 
that of bulk systems (for small d) and approaching those 
of monolayers (for large tube sizes). 

The electronic structures of M0S2 and WS2 and the 
resulting optical properties come from the d-electron or- 
bitals that dominate the valence and conduction bands 
(see Fig. [5] for M0S2). The projected DOS (PDOS) of 
TMS2 shows that p-states of sulfur atoms hybridize with 
the d-states of the transition metal atoms at the top of va- 
lence band and the bottom of conduction band. The core 
states are dominated by the s-orbitals of the chalcogenide 
atom. The PDOS of tubular structures do not change sig- 
nificantly from that of layered forms. 



4 Conclusions 



We have studied electronic properties of TMS2 (TM 
— Mo and W) in the layered and tubular forms. For the 
layered structures we have found the crossover from an 
indirect- to a direct-band gap semiconductor, when going 
from bulk to a monolayer limit. The strain energies, orig- 
inating from rolling up a monolayer to a tube, are one 
order of magnitude larger than in case of SWCNTs, but 
decrease for larger tube diameters. These strain energies 
scale as 1/d^ {d - tube diameter), confirming also the re- 
suhs of Seifert et &\.^ obtained at the DFTB level. 

There are no size dependent electronic irregularities 
of the studied TMDC-NTs, as it is found for SWCNTs. 
The band gap values increase with increasing the tube 
size and range from that of a bulk structure to that of a 
monolayer. Zigzag NTs have larger band gaps than arm- 
chair ones for a given d. The very interesting phenomenon 
was found for the electronic band structures: zigzag NTs 
resemble the band structure of monolayers, while the arm- 
chair NTs (especially for large tube sizes) resemble the 
electronic structure of bulk systems. 
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Fig. 6. (Color online) Band structures of tubular TMS2 (TM 
= Mo and W) with armchair and zigzag arrangement and chiral 
index n = 11, 15 and 24. The arrows indicate the fundamental 
band gap (direct or indirect) for a given system. The top of 
valence band (blue) and bottom of conduction band (green) 
are highlighted. 
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Fig. T. (Color online) Calculated gap energies (Zi) of (n,0) and 
(n,n) TMS2-NTS (TM = Mo and W) with respect to the tube 
diameter (d). The corresponding layered structures are given 
as horizontal lines: solid-bulk, dashed-monolayer. 



MoS^ monolayer 




Eo [eV] 



Fig. 8. (Color online) Partial density of states of bulk and 
monolayered M0S2. The projections of Mo and S atoms are 
given together with the contributions from 4d and 3p orbitals 
of Mo and S, respectively. The vertical dashed lines indicate 
the the top of valence band. 
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